Atoms in a radiofrequency-dressed optical lattice 
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We load cold atoms into an optical lattice dramatically reshaped by radiofrequency (rf) coupling 
of state-dependent lattice potentials. This rf dressing changes the unit cell of the lattice at a 
subwavelength scale, such that its curvature and topology departs strongly from that of a simple 
sinusoidal lattice potential. Radiofrequency dressing has previously been performed at length scales 
from mm to tens of /im, but not at the single-optical- wavelength scale. At this length scale significant 
coupling between adiabatic potentials leads to nonadiabatic transitions, which we measure as a 
function of lattice depth and dressing frequency and amplitude. We also investigate the dressing by 
measuring changes in the momentum distribution of the dressed states. 

PACS numbers: 03.75.Mn, 03.75.Lm 



For almost twenty years, optical lattices have been 
used to cool and confine neutral atoms, recently lead- 
ing to condensed-matter-analog systems, optical clocks 
and possible platforms for quantum computing |2 12 
El m |5l E]. The appeal of optical lattices arises from 
their versatility and tunability, e.g., through control of 
the topology, depth, and spin dependence of the lat- 
tice [II m [H [To]. Here we extend this control by rf 
dressing a spin- dependent lattice, a technique which has 
been recently proposed [T^. The dressed lattice com- 
prises unit cells whose localized adiabatic eigenstates are 
spatially- varying superpositions of the bare spin states. 
The adiabatic potentials have new subwavelength struc- 
ture, which in certain limits can be ring-like within a unit 
cell. 

The creation of rf-dressed adiabatic potentials in inho- 
mogeneous magnetic fields is well-known experimentally, 
from the use of rf transitions in evaporative cooling [12^ 
to the creation of shell-like quasi-two-dimensional sys- 
tems in magnetically-trapped condensates [T3l [TH |T5], 
and the use of rf dressing to effect condensate splitting 
and interferometry in micron-scale chip traps [16 . Op- 
tical lattices with sub-half-wavelength structure in ID 
have been generated by several groups through the use 
of multiphoton Raman processes [171 [HI [E] oi" static 
magnetic field couplings [20 . The approach we present 
here, using the rf dressing of a state-dependent lattice, 
creates 2D subwavelength structure beyond simple lat- 
tice period division. We create such a dressed lattice, 
and load ultracold atoms into the ground band of the 
uppermost adiabatic potential. 

For these state-dependent lattice experiments we use 
a spin-polarized ^^Rb (5^ 51/2 |^ = I^ttif = — 1)) Bose- 
Einstein condensate (BEG). The state dependence is pro- 
duced by tuning the lattice laser between the ^^Rb Dl 
and D2 resonances, where vector light shifts are im- 
portant [9 . Applying an rf field of appropriate ampli- 
tude and frequency creates an adiabatic potential whose 
shape is distorted from the bare potentials, as shown in 



Fig. [TJb). After holding the atoms in the dressed lat- 
tice for a variable time, we abruptly turn off the lattice 
and observe the momentum distribution by imaging the 
atoms after time- of- flight (TOF) time txoF- This gives 
information about the structure of the dressed lattice, as 
well as a measure of the lifetime of the dressed adiabatic 
eigenstates. 

Our apparatus, described elsewhere [ST, produces 
BECs of ~10^ atoms in a loffe-Pritchard magnetic trap. 
We load the = — 1 condensate into the ground band 




FIG. 1: (color online) (a) The energy levels of ^^Rb relevant to 
our lattice; the rf coupling can selectively address either mp 
transition due to the quadratic Zeeman shift, (b) Top two 
images: the bare lattice potentials for mp = ±1, exhibiting 
approximately opposite light shifts. Bottom two images: the 
uppermost adiabatic potential, weakly dressed and strongly 
dressed, (c) From top: bare lattice potentials for U ^ IOEr; 
dressed adiabatic potentials at a;rf/27r = 35.90 MHz, Q/27v = 
20 kHz; dressed adiabatic potentials at a;rf/27r = 35.90 MHz, 
Q/27T = 205 kHz. Note the significantly increased gap for 
stronger coupling. 
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of a 3D optical lattice via an exponentially increasing in- 
tensity ramp of duration 300 /is (r = 50 /is) , a timescale 
adiabatic with respect to band (vibrational) excitation. 
We then remove the magnetic trapping fields, leaving a 
uniform bias field of 5.117(3) mT [27^, a Zeeman res- 
onance = 36.12(2) MHz and a quadratic Zeeman 
shift b' j2i\ = i^_i^o — ^o,+i = 376(1) kHz, where hvrn.m' is 
the positive energy difference between Zeeman sublevels 
\F = l,mF). 

The lattice in the x-y plane {xy lattice) is similar to 
that described in [21 . We tune a Ti:sapphire laser to 
A =790.76 nm, where for cr^ light, atoms in the three 
F = 1 Zeeman sublevels tuf =— 1, 0, +1 experience 
total light shifts in the approximate ratio of —3 : 1 
: 5. This can be understood as the sum of a state- 
independent scalar light shift and a state-dependent vec- 
tor light shift — an effective magnetic field related to the 
local ellipticity of the optical polarization. For appropri- 
ate phase shifts of the lattice beams [9 , the combination 
of the light shifts and the bias field along {x — y) / V2 
results in the potentials shown in Fig. [l] 

The atoms are confined along z (the direction of grav- 
ity) using a ID optical lattice, derived from the same 
Ti: sapphire laser [28]. These lattice beams are lin- 
early polarized, yielding an effectively blue-detuned and 
spin-independent lattice. As measured by pulsed-lattice 
diffraction [22] this vertical lattice is QEr deep, where 
Er h^k'^/2M ^hx 3.68 kHz, k 2ti/\ and M is the 
atomic mass. This lattice, together with the xy lattice, 
fully confines the atoms in 3D, and we do not observe 
any displacement due to gravity over the duration of the 
experiment, even though the magnetic trap has been re- 
moved. 

After loading, we wait 30 ms for the atoms to dephase 
and fill the lowest Bloch band of the lattice [8 , such that 
a sudden-release TOF image will approximate the single- 
lattice-site momentum distribution (with resolution lim- 
ited mainly by the spatial extent of the atom cloud). The 
dephasing time is well within the 1/e lifetime of atoms 
in the lattice due to spin-fiip loss from light scattering, 
which we measured to be 85 ms at an xy lattice depth 
oi U = 55Eji for ttif = — 1 atoms; this lattice depth is 
more than twice the largest used in the experiment. The 
depth of the xy lattice was measured with two-photon 
Raman vibrational spectroscopy [23]. 

We dress the lattice with an rf magnetic field aligned 
perpendicular to the bias field. The single-particle 
rotating- frame Hamiltonian, taking the rotating-wave 
approximation (RWA), is H = p^/2M + Hi(r) , where 
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and S = LOrf — 27ri^_i^o • The state-dependent lattice 
potentials Vmpi'^) (see Fig. 1) are calculated using our 



FIG. 2: (color online) (a) The width of the momentum dis- 
tribution as a function of final dressing frequency, as mea- 
sured by absorption imaging. The width is calculated as 
Md/ Mtof, where d is the Gaussian 1/e radius along one 
direction, corrected for the initial cloud size. The displayed 
uncertainties represent la shot-to-shot scatter of five points. 
The shaded area represents predicted widths given the values 
and estimated uncertainties of the magnetic field (5.117(3) 
mT) and the depth of the bare lattice U = 10.0(5)^^^. Also 
shown are the bare and dressed single-site potentials (solid 
line) and ground-state wavefunctions (gray). (b,c) Sample 
images of the off-resonant (effectively bare) and dressed mo- 
mentum distributions respectively (both at Q/27r = 205 kHz) 
showing significant narrowing of the latter. (d,e) Examples of 
dressed momentum distributions created through faster (> 1 
MHz/ms) rf ramps, such that higher vibrational states of the 
dressed lattice are populated. 



model of the lattice [21]. Q is obtained through obser- 
vations of rf-driven oscillations between tuf = —I and 
= (in a spin-independent lattice) at rates up to a 
maximum of cJosc/^tt = 200 kHz. Q differs slightly from 
cjosc due to the multilevel nature of this system. Diago- 
nalizing only Hi (i.e., the Born-Oppenheimer approxima- 
tion) gives rise to adiabatic potentials; failure of this ap- 
proximation will manifest as momentum-dependent cou- 
plings between the adiabatic potentials. 

We apply rf beginning well below resonance (::^ — 1 
MHz), sweeping at constant rate (300 kHz/ms) over :^ 
3 ms to a final near-resonant frequency, and hold the 
newly created dressed state for a variable time. The lat- 
tice beams, rf field, and bias field are then turned off in < 
1 /is, 10 /is, and 300 /is, respectively. After TOF of 
12.2 ms we observe the atoms with resonant absorption 
imaging along the z direction. 

Fig. [2ja) depicts the observed width of the momentum 
distribution (from a Gaussian fit to the TOF density dis- 



3 



tribution) as a function of final frequency at maximum 
rf power. The optical power in the xy lattice was fixed 
such that the bare depth U was :^ I^Er. Also shown 
is the expected width from the calculated ground state 
momentum wavefunction in the adiabatic potential. The 
data show narrowing of the momentum distribution as 
the lattice is dressed, implying significant alteration of 
the wavefunction and hence the lattice structure itself. 
At these dressed depths, while the potential exhibits ring- 
like character, the ground state is still simply connected: 
see Fig. 2(a). The observed dressed state (corresponding 
to rriF = —1 without rf) is predominantly a superposi- 
tion of rriF = ±1 at these rf detunings [29 . Unlike in a 
two-level system, the avoided-crossing energy gap HA be- 
tween adiabatic potentials is not simply related to Q. For 

< ^^ A n'^/S\ and for Q ^ , A n/V^ - d'/4.. 
For Q/27T = 205 kHz, the largest used in this experiment, 
and c^rf = 35.90 MHz, we calculate A/27r :^ 75 kHz, a 
value only weakly dependent onU. To further character- 
ize our dressed lattice, we performed Raman vibrational 
spectroscopy, but observed that the Raman signal did 
not persist beyond weak dressing. This may be due to 
to shrinking matrix elements between the combined spin 
and spatial eigenfunctions of the dressed bands, or in- 
creasing width of the excited bands of the (shallower) 
dressed lattice. 

Three adiabaticity criteria are relevant to loading and 
holding atoms in the ground band of the dressed lattice: 

(1) Adiabatic following of the local spin eigenstates of 
7^1 (r) during the upward rf sweep S{t) is well-satisfied. 

(2) Deformation of the dressed lattice with respect to 
vibrational excitation is nearly adiabatic, yielding pre- 
dominantly ground-band occupation. Faster sweep rates 
yielded interesting deviations from this condition, shown 
in Fig.|2|d,e), in which we observe lobed structure in the 
momentum distributions indicating higher vibrational 
states. (3) The Born-Oppenheimer approximation, i.e., 
the degree to which zero-point motion in the lattice does 
not induce transitions between adiabatic potentials, is 
more difficult to satisfy. In previous experiments involv- 
ing rf-dressed adiabatic potentials, the length scales were 
large enough such that this condition was easily met for 
cold atoms [131 [Hi [151 |16J . For our experiment, the bare 
lattice confinement is such that, even for our largest rf 
coupling, nonadiabatic loss is a factor; in our experiment 
these loss rates are 10^-10^ s~^. These losses do play a 
role in the loading process itself, but this is irrelevant as 
long as enough atoms are loaded such that subsequent 
loss is measurable. 

We load atoms into the uppermost adiabatic dressed 
potential using the procedure described above, and mea- 
sure the number of atoms remaining in the central mo- 
mentum feature, as in Fig. 2(b-e), as a function of hold 
time. Atoms transitioning to lower adiabatic potentials 
appear in absorption images as high-momentum rings of 
~ lOhk (well outside the range of Fig. 2(b-e)), and are 
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FIG. 3: (color online) Dressed- state losses as a function of 
calculated avoided-crossing gap A/27r, for four bare lattice 
depths. Vertical bars represent uncertainties-of-fit; multiple 
points near a given A indicate scatter in multiple runs. The 
data are well-represented by 7(A) = Ae~'^^^'^^ , with B — 
83(6) /is , 83(11) /is, 82(13) /is , and 100(30) /xs, for U = SEr, 
IOEr, 12Er, and 16Er respectively. The scaling of the data 
does not change if plotted versus Q. 



counted as loss. Fig. [3] shows loss rates as a function 
of calculated avoided-crossing gap A as a function of /7, 
at dressing frequencies near 35.90 MHz. The rates are 
first obtained through two-parameter fits to an exponen- 
tial decay and then slightly corrected via subtraction of 
the appropriately scaled spin-fiip loss rate (1/85 ms~^ for 
55Er). The data are well-described by 7(A) = Ae~^^. 
Figure |4] shows loss rate at constant as a function of 
/7, wen described by ^{U) = Ce^^ . 

We consider our observations in the context of semi- 
classical theory for the losses. The traversal of a two- 
level avoided crossing is often described by Landau- 
Zener (L-Z) theory, which for a particle traveling along 
X at velocity v through an avoided crossing character- 
ized by a gap A and bare energy levels whose differ- 
ence E{x) = Ex^ the adiabatic following probability is 
Pa = 1 — e~h^^^ /(^^ ). A semiclassical scaling argu- 
ment for the motion of an atom in the ground state of a 
lattice site of depth U yields v cx /7^/^ and vE' oc /7^/^, 
and a loss rate scaling as 7 cx cc;^(l — Pa), where uoi is 
a trapping frequency and (1 — Pa) <C 1. For a parti- 
cle of insufficiently large energy with respect to A (i.e., 
atoms trapped in the ground state of an adiabatic poten- 
tial), the simple L-Z approach may not apply. In a more 
sophisticated semiclassical approximation 125] , more 
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FIG. 4: (color online) Dressed-state losses as a function of 
bare lattice depth [/, at constant Q./2n = 205 kHz and dress- 
ing frequency 35.900(35.875) MHz for the first five (last four) 
depths, resulting in approximately constant gap A/27r 75 
kHz. Vertical bars represent uncertainties-of-fit; multiple 
points at a given depth represent the scatter in multiple runs. 
Horizontal bars represent known drifts in laser power; not 
shown is a constant systematic uncertainty in depth determi- 
nation on the order of 5%. The data are well-represented by a 
simple exponential fit ^{U) = Ce^^, yielding C = 1.2(4)s~^ 
and D = 0.27(3)^^\ 



ficient rf coupling; based on the data presented here, we 
extrapolate that ringlike wavefunctions of lifetime :^ 100 
ms should appear for a coupling Vt/2T: 400 kHz and 
U 55 Er. 

This paper presents observations of a novel lattice po- 
tential consisting of rf-coupled components of a state- 
dependent optical lattice. The lifetime of this dressed 
lattice depends exponentially on both coupling strength 
and optical lattice depth. We observe lifetimes up to 100 
ms for strongly-dressed shallow lattices, a timescale suffi- 
cient for the realization of many-body condensed-matter 
systems. Extensions of the band-structure engineering 
demonstrated by this experiment lead to lattice geome- 
tries with structures much more complex than possible 
with conventional optical lattices, leading to intriguing 
single-site wavefunctions and the next-nearest-neighbor 
interactions required for extended Bose-Hubbard mod- 
els [26 . 

We thank P. Zoller for inspiring this particular exper- 
iment and A. J. Daley, G. Pupillo, and P. Z. for ongoing 
helpful discussions. This work was partially supported 
by DTO and ONR. N.L., P.J.L., and B.L.B. acknowl- 
edge support from NRG fellowships. 



appropriate to ground-state trapped atoms, the result- 
ing transition rates can be estimated as 7 oc a;^e~"^/'^% 
where ui oc y^U/A is the trapping frequency {ui <C A) 
of the uppermost adiabatic potential, and a ~ 1 is a 
constant. The shape of our adiabatic potentials differs 
significantly from this simple case, but the scaling sug- 
gests transition rates that change exponentially with the 
coupling gap and bare lattice depth. 

Unsurprisingly, the L-Z model does not describe the 
data in Figs. |3] and |4j the exponential coefficient scales as 
—A rather than — A^ and also scales as [/, not 
Fits to the L-Z model are significantly worse than those 
in Figs. 3 and 4, although they agree with the general 
trends. Similar results were found upon calculating a 
Fourier- weighted sum of the L-Z formula over the dressed 
2D ground-state wavefunction, J |?/^kP^NA(k)(ik. The 
more sophisticated semiclassical approximation performs 
better than the L-Z approach, but still fails to satisfac- 
torily predict the observed scaling. It is likely that ID 
semiclassical descriptions are simply inadequate for this 
2D quantum problem; perhaps a more complete solution 
will yield the very simple scaling we observe. Further- 
more, heating to higher vibrational levels of the upper- 
most adiabatic potential, likely caused by instability in 
the rf dressing frequency (measured to be on the order 
of several kHz), is likely to increase the nonadiabatic loss 
rate due to increased mean velocity in higher bands. 

While the dressed potentials shown in Fig. [TJb) are 
ringlike, the corresponding ground-state wavefunctions 
are generally not. Nevertheless, ringlike wavefunctions 
with reasonable lifetimes would be possible had we suf- 
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